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Abstract 

We discuss the function Wt(x) defined via the implicit equation Wt(x)-tan [W t (x)] = 
x, which appears in certain quantum mechanical and field theoretic applications. 
We investigate its analytic structure, develop series expansions for both small and 
large x, and provide various techniques for its numerical evaluation in the complex 
plane. 
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1 Introduction 



Implicit equations for functions appear frequently in physics although they are not dis- 
cussed at all in standard handbooks on (standard) special functions [1, 2]. The most 
famous example, probably, is Kepler's equation 

out = ip — e sin-?/' (1.1) 

which describes the dependence of the angle ip ( "the eccentric anomaly" ) on the time t 
for an elliptical orbit with eccentricity e [3] . It has been a challenge to solve this equation 
for generations of astronomers and mathematicians, among them F. W. Bessel who was 
led by this problem to the functions which bear his name. 

Another transcendental equation which has made its appearance in a variety of appli- 
cations is 

W(x) exp[W(x)] = x, (1.2) 

the history of which goes back to J. H. Lambert and L. Euler in the 18th century. A 
comprehensive survey of the many available results for this function has been given in 
Ref. [4] where it has been called the Lambert W-function (probably because Euler already 
has enough named relations and functions to his credit Apart from the examples 

given in the above article there is also some recent interest in the properties of W(x) in 
particle [6] and nuclear physics [7]. 

Here we want to study a function defined by a superficially similar equation as given 
in (1.2), viz. 

W t (x) • tan [W t (x)] = x (1.3) 

which, however, has a much richer structure than Lambert's W-function. Because of the 
similarity we propose to call 2 the solution of Eq. (1.3) W t (x). Analogously, one could 
define 

W rt (a) -cot[W ct (x)] = x (1.4) 

and W s (x) and W c (x). Eq. (1.3) emerged in a variational approach to Quantum Elec- 
trodynamics [8] where the relation between the bare mass and the physical mass of the 
electron was obtained in terms of the solution of an implicit equation of the type (1.3). 
However, as we will see in Section 2, the same equation also appears as an eigenvalue 
equation in a one-dimensional quantum mechanical problem for an infinite square well 
with a residual (^-function interaction. If the strength of this additional interaction is 

1 In Ref. [5] a slightly more general equation is called "Ramanujan's equation". 

2 "Names are important" [4] for a consistent nomenclature. In the present scheme the Lambert W- 
function would be denoted by W e (x). 



2 



constant one is naturally led to Eq. (1.4) whereas a linear energy dependence of this 
strength leads to Eq. (1.3). 

In the present, mostly elementary, note we will concentrate on the W t -function and, 
since we are not aware of any other systematic study, we will derive some representations 
by simple methods and study the interesting analytical structure of this function in the 
complex plane. We hope that this collection of results will be helpful for other applications 
where this function makes its appearance. 



2 A quantum mechanical problem leading to W t (x) 

Consider a nonrelativistic particle moving in an one-dimensional infinite square well in 
< £ < a. Assume that an attractive 5-function interaction is added in the middle of the 
interval whose strength is growing with the energy of the particle 3 

AV = -\E6(Z-?) . (2.1) 

As the total Hamiltonian is symmetric with respect to transformations £ — > a — £ the 
wave functions can be classified according to their symmetry under this transformation 
^n(0 — =tVVi(a ~ £)■ Only the even wave functions are influenced by the additional 
interaction (2.1) since the odd wave functions vanish at £ = a/2. The Schrodinger equation 

2m d£ 2 

with E = h 2 k 2 /(2m) has the solutions 



XEs((-^j m = E-m (2.2) 



A n sm(k(a — £)) , a/2 < £ < a . 

The solution (2.3) fulfills the boundary conditions that the wave function has to vanish 
at the edges (£ = and £ = a ) of the infinite square well. At £ = a/2 the wave function 
has to be continuous 

A/Sin (y) = A//Sin (y) ' (2 - 4) 

but its derivative makes a jump proportional to the strength of the 5-function. This can 
be derived in the standard way by integrating Eq. (2.2) from a/2 — e to a/2 + e and letting 

e -> 

~ r dip ( a \ dip f a 



2m 



di \2 + V di \2 



-AA'r(|) = 0. (2.-T) 



3 Energy-dependent potentials are frequently considered in nuclear and hadronic physics where they 
arise from eliminating part of the physical Hilbert space. 
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The r.h.s. vanishes since the wave function is finite and therefore its integral over an 
infinitesimal distance does not give a contribution. For our solutions this implies 

( ka\ t , ( ka\ 2m, ^ t ( ka\ , 9, . ( ka\ 
A u k cos I — J +Ajk cos I — j = j^X E A 7 sin I — J = A A/ sm I — j . (2.6) 

Combining the boundary conditions (2.4) and (2.6) we obtain as eigenvalue condition 4 

ka ( ka\ a , „. 

with the general solution 

* = fw,(i). (2.8) 



X 



-3tc/2 



i-tu/2 



-2- 



-4- 



2 



4 



Figure 1: Eigenvalue equation (2.7) with x — a/X and y = 



The graphical representation of the eigenvalue equation in Fig. 1 shows that W t is a 
multivalued function, and we have to define its branches. Our definition will be motivated 
by physical arguments: for vanishing interaction (A = 0) the r.h.s. of Eq. (2.7) diverges 
and thus the tan-function has to diverge as well. Therefore 

ka 

T 



A=0 



= ±- ±- 

2' 2 ' 



(2.9) 



4 Note that the energy-dependence of the additional interaction is essential; without that one would 
obtain the relation (ka/2) cot(fca/2) = \am/(2h 2 ) leading to the W c t-function. 
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which obviously gives the unperturbed wavenumbers (and hence energies) of the symmet- 
ric states in the infinitely deep square well of length a. Thus we define the nth branch of 
Wt by the requirement that x = oo is a regular point and 



W[ n \x — > oo) = sgn(n) — ^ ir , n = ±1, ±2, ... 



(2.10) 



The limit of very strong coupling (A — > oo) is also easily understood: if the (^-function is 
attractive the energy levels are lowered: 



ka 

T 



A^+oo 



0, ±7T, 



(2.11) 



i.e. 



W[ n \x -> +0) = sgn(n) (\n\ - 1) tt , n = ±1, ±2, . . . . (2.12) 

However, if the 5-function is very much repulsive the square well is practically halved and 
the eigenvalue condition is sin(ka/2) = 0. This is formally the same as for A — > +oo but 
k = is excluded as trivial eigenvalue. Therefore the levels are pushed up: 



ka 

Y 



A^ — oo 



±7T, ±27T, 



(2.13) 



i.e. 



W^ n) (a: 



-0) = mx . (2.14) 

These conventions are summarized in Fig. 2. Since the value of Wt(0) depends on the 
way one approaches the point x — 0, it must be a branchpoint for the function. Indeed, 
replacing for small argument the tan-function by its argument one finds for the principal 
branch (which corresponds to the ground state of the system) 



W t (1) (x) 

In general we have from eqs. (2.12, 2.14) 



y/x . 



w. 



(n), 



-o) 



W t (n+sgn(n)) (+0) 



(2.15) 



(2.16) 



and our branch conventions can be summarized by writing the defining Eq. (1.3) in the 
form 



W[ n \x) = sgn(n) (jn| - ^ tt + ^ In 



x — iW[ n ^ (x) 



x + iW 



(»»), 



X 



sgn(n) ( |n| — ^ ) tx + &(—x) sgn W^(x) 



(2.17) 

x-iW{ n \x)p.l8) 



ii + arg 

where ln(z) denotes the principal branch of the logarithmic function (— n < arg(^) < tt). 
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n=+1 
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-jc/2 
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Figure 2: Branches of the W t -function. 



From Eq. (2.18) we obtain the important relation 



W 



(-«), 



2 

= — sgn(n) (\n\ — - ) n — 0(— x) sgn 

= -M n \x) 



(x) = — sgn(n) ( \n\ — ^ ) n + 0(— sgn w| ^(x) 7r + arg 



x — iW 



(-n) 



-w t ( - n) (x) 



7r — arg 



(2.19) 



which allows us to restrict our analysis to branches with n > (cf. Fig. 2). It is, of 
course, easy to find some special values of the W t -function, e.g. 



Wt (n+1) Q + wtt) = ^ + nvr, n = 0, 1, ... 



(2.20) 



The above definition of the branches of w[ n \x) is still incomplete because we have not 
yet specified how different sheets of the Riemann surface are connected to each other 
across the cuts that must also be explicitly defined. This will be discussed in more detail 
in Section 5. In the next two Sections, our discussion concerns only the properties of the 
main branch (x) on the real x axis, where (x) is continuous except at the point 
x = 0. 
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3 Series expansions, differential equation and inte- 
grals 

In this Section we concentrate on the the principal branch of the W t -function, i.e. n—1. 
From the defining equation (1.3) it is easy to obtain the first terms in the series expansion 
for small x > 



W^(X) = = 



fc=0 



1 1 11 

1 x H 

6 360 



x 



17 
5040 ~ 



281 



-x 



604800 



x* + 



(3.1) 



and, by setting W t = vr(l — v)/2 for large 
W«(x) 



— OO L 

z fc=0 x 



7T 

2 



7T 



X X' 



12 x 3 



- 1 



+ 



X 1 



(3.2) 



To calculate the expansion coefficients for \x\ — > oo systematically, one may use Lagrange's 
expansion theorem for implicit functions 5 : 



b = 1 , b k = 



1 d 



k-l 



v(l — v) 
2 J k\ dv k - x \ v tan(7rt;/2) / 



, k > 1 . 



(3.3) 



d=0 



The convergence radii of the above two expansions may be estimated from the root test 



x < Po 



lim | cfcfc j 



lim Pq 



(k) 



\x\ > p c 



lim \bk 



1/fe 



lim p« . (3.4) 

fe^oo 



Fig. 3 shows the result from the first 100 coefficients; it seems that these expansions 
converge only for finite values of x < 2.8 and l/\x\ < 2.4 which hints at additional singular 
points in the complex x-plane. We will locate these branchpoints of the W t -function in 
Section 5. 

At high orders in the above expansions the computational cost of using Lagrange's 
expansion soon becomes prohibitive. Indeed, the data in Fig. 3 were generated by develop- 
ing recursion relations for the expansion coefficients in Eqs. (3.1, 3.2). By differentiating 
Eq. (1.3) with respect to x and eliminating the trigonometric functions one arrives at the 
simple differential equation 



W t '(:r) 



W t (x 



(3.5) 



x + x 2 + W t 2 (x) ' 

Turning first to the expansion around \x\ — > oo, we set t = 1/x and rewrite the above 
differential equation as 



d ( 1 \ „ dW t (l/t) 
[ + } dt \W t (l/t)) dt 



W t (l/f) 



(3.6) 



5 See, e.g., Ref. [1], Eq. 3.6.5 . 
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Figure 3: Estimates of the convergence radius for the series expansions (3.1, 3.2) oiW[ 1 \x) 
from the first 100 coefficients. 



With the ansatz 

one obtains a system of coupled recursion relations for the expansion coefficients b k , c k 

(k + l)c k+1 + (k-l)c k = ^(1-<W (k-l)b k -! (3.8) 
k 

cj = 5 k0 , (3.9) 

j=0 

which can be solved successively for k = 0, 1, . . . starting with bo — 1. For asymptoti- 
cally large values of k, we find (numerically) that the coefficients b k approach the simple 
functional form 

bk k ~^* ~j^j2 (P°°) k sin ( ak + b ) > ( 3 - 10 ) 

with c ~ 0.584, a ~ 2.25, b ~ —3.59 and the radius of convergence of the expansion in 

t is 1/poo — 1/2.64 6 . Indeed, by writing the trigonometric function as a sum of two 

6 In the quantum electrodynamical context where the W t -function appeared [8], this implies a finite 
convergence radius of the perturbative expansion contrary to the usual semi-heuristic arguments [9] for 
an essential singularity at zero coupling constant. 
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exponentials we find that we can absorb the rapid oscillations in b k into the definition of 
a new variable x' (or, equivalently, t'), obtained by rotating x to xe ta (or xe~ ta for the 
second of the two exponentials). Hence we anticipate that the branchpoints responsible 
for the finite radius of convergence of W t are in fact located at x — Pooe ±m ~ 2.64e ±a727 ", 
i.e. in the second and third quadrants of the complex x plane. Furthermore note that, 
generically, the expansion coefficients of the root (1 — tp^Y behave, at large orders k, 
like (poo)^ /k K+1 , which leads us to suspect that the branch cuts are of square root type. 
These suspicions will be confirmed in Section 5. 

For x — > +0 on the other hand we set w[ 1 \x) = y/xf(x) and obtain the differential 
equation 

(x - l)f(x) + 2x(x + l)f(x) + f(x) + 2 -x^ (f(x)) = . (3.11) 

The ansatz 

oo 

f(x) = Y, d kX h (3.12) 

fc=0 

now leads to the following system of recursion relations for the coefficients 

1 2k + 3 

a k + (1 - 5 k0 ) a fc _i + - 2fe _ d k = (3.13) 
1 ^ 

- ^ (4j - k) a, d k -j = 4, k = l,2,... (3.14) 

where ao = do — 1 is the starting point. Eq. (3.14) follows from the identity / • (/ 3 )' = 
3/ 3 • /'. For large k the recurrence relations yield 

ak ~j^j2 (Po)' h sin(ak + b) , (3.15) 

with c ~ 0.564, b ~ —3.14, p — Poo and a the same as above. 

Finally, we note that many integrals over the W t -function may be performed analyti- 
cally by a partial integration and changing variables to x — ytany: 

Jdxf[W t (x)} = xf[W t (x)} - Jdyytfmyf'iy) . (3.16) 



For example, 



/ 



dx\n[W t (x)} = x\n[W t (x)] + In 



/ 



cos(W t (a;)) 
1 



(3.17) 



dxln[sin(W t (a;))] = a: In [sin (W t (x)) ] - -W t >) , (3.18) 



and so on. 

From the latter equation and the asymptotic behaviour of the W t -function one immedi- 
ately obtains the definite integral 

r 2 



I 



dx In 

o 



sin(W t (1) (x))l = . (3.19) 
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In addition one finds 

/""V 4 mi 7T 2 7T 1 

/ drW 1 '^) = - + -In2--G (3.20) 
Jo lb 8 2 

where G = 0.91596594... is Catalan's constant. This is just the special case when Young's 
inequality (Ref. [2], Eq. 12.315) applied to the integral 

f dxW^ix) > a- + -\n2- -G , a< - (3.21) 
Jo tw ~48 2 ' - 4 v; 

becomes an equality due to the relation (2.20). 



4 Numerical evaluation and approximations to the 
W t (x)-function 

4.1 Chebyshev approximation 

The numerical evaluation of the W t -function for real arguments is most easily accom- 
plished by expanding in Chebyshev polynomials [10] 

[ v^E fe =o« fe T fe (f -l) , 0<x<a 
W t (1) (x) = f E k =oPkT k h) , |x| > a (4.1) 

[ vrE fc =o7feT fc (f + l) , -a<x<0 . 

Table 1 gives the corresponding coefficients for a = 3.5, which was found to be a convenient 
value to use in order to divide the different regions. One should keep in mind that 
|^fi(^)| < 1 for x £ [—1,-1-1], so that (assuming a rapid decrease of the coefficients) the 
absolute value of the last coefficient gives an estimate of the accuracy with which (x) 
has been approximated by the truncated Chebyshev expansion. 

4.2 Approximation via iteration 

To calculate the W t -function for arbitrary (complex) x, one may naturally try to solve 
Eq. (1.3) by iteration. In this respect we have found Halley's improvement [11] to the 
Newton-Raphson algorithm to be useful. After % iterations, this leads to the estimate for 
the root y of a function f(y) as 

(4.2) 



Vi+i = Vi - f{Vi) 



f'iVi) 



f(vi) r (y» 
V'ivi) 
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k 




Pk 


7fc 





0.80600536 


1.03465858 


0.82312766 


1 


-0.16766125 


-0.28291110 


0.16771494 


2 


0.02302848 


0.03258714 


0.01423939 


3 


-0.00298934 


0.00177957 


-0.00442520 


4 


0.00030980 


-0.00206359 


-0.00095545 


5 


-0.00001275 


0.00044900 


0.00024196 


6 


-0.00000478 


0.00002021 


0.00007926 


7 


0.00000178 


-0.00004188 


-0.00001706 


8 


-0.00000038 


0.00001127 


-0.00000743 


9 


0.00000006 


0.00000018 


0.00000136 


10 


-0.00000001 


-0.00000111 


0.00000075 


11 


0.00000000 


0.00000035 


-0.00000012 


12 


0.00000000 


-0.00000001 


-0.00000008 


13 


0.00000000 


-0.00000003 


0.00000001 


14 


0.00000000 


0.00000001 


0.00000001 



Table 1: Chebyshev expansion coefficients of Wj; '(x) in Eq. (4.1) for a = 3.5. 



which, with f(y) = x — yt&ny and x fixed, yields 



Vi+i = Ui+(x- Vit&nyi) 



yi(l + tan 2 y { ) + tan^ 



(Vitanyi + l)(tan 2 ^ + 1) n ~ 

yi(l + tan i/J + tanyj 



(4.3) 



4.3 A variational principle for W t (x) 

The quantum mechanical model of Section 2 can be used to obtain some insight and some 
approximations for the W t -function. It is slightly nonstandard because of the energy- 
dependence of the additional interaction but this energy-dependence is linear so that it 
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can be easily handled. The Schrodinger equation may be cast into the form 



ffolKO = EN(0m (4.4) 

with Hq = —h 2 d 2 /(2md^ 2 ) and iV(£) = 1 + A5(£ — a/2). For a positive definite operator 
N we can define 

H = N-^HvN- 1 ' 2 , = N 1 ' 2 ^ (4.5) 

and end up again with a normal Schrodinger equation H<f) = E<f> for which all standard 
methods of quantum mechanics can be applied. In view of the finite convergence radius 
of the series expansion the variational method is of particular interest here as it is a 
non-perturbative method. It gives an upper bound for the true energy 

E= ™ < <^l«k>, (4.6) 

where (f> t is a trial state containing variational parameters. Using Eq. (4.5) this can be 
cast into the form 

E < <^ g °l^> > (4.7) 
- < Wt> 

which is a well-known variational principle for the generalized eigenvalue problem (4.4) 
[12]. Recalling Eq. (2.8) and setting a = n for convenience we obtain 



/■7T / //, i / ( \ \ ,-r: — 7T 



,^(0)=^W = 0. (4i 



As usual this variational principle can be applied most easily to the ground state, i.e. to 
the principle branch. Indeed, taking as trial wave function just the unperturbed ground 
state wave function ^(0 = srn £ gives 



which covers both limiting cases of the series expansion studied in the previous Section 7 . 
However, it predicts a square-root type branchpoint on the real axis at x = — 2 (where 
< ip t \N\ip t > vanishes) which obviously is not realistic. One can easily do better by 
allowing for an admixture of the first excited symmetric state, i.e. ipt(Q = sin£ + b sin 3£, 
where b may be optimized for each value of x. A simple calculation in which we fix the 
square root sign by demanding that b — > for x — > ±oo then gives 

m i] (x) < —( y /2 . (4.io) 

2 \5x + 10 + 2 sgn(x) V25 + 16a; + Ax 2 ) 
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^ 4 

-t— ' 

3 



2 

















*\ variational 




exact 


% 

\ I 
\i 
\l 


— 


1 1 1 1 — 





-4 -2 2 4 



1/X 

Figure 4: The principal branch of the W t -function as function of l/x, as well as its 
variational approximation of Eq. (4.10). 

This improves the bound near x = +0 to W^ix) < 3nV^/20 = 1.0537^/5 and predicts 
W; i} (-0) < ttv/5/2 = 1.1180 7T. 

Fig. 4 compares this variational approximation with the exact Wj^-function calculated 
numerically via the iteration procedure. For x > the agreement is excellent. While the 
approximation is regular on the negative x-axis it exhibits branchpoints at x — — 2±3i/2 
in the complex x-plane where the inner square root vanishes. This would restrict the series 
expansions of (x) around x = or x = oo to converge for magnitudes smaller or larger 
than 5/2, respectively. In the next Section we will find that the exact convergence radius 
is 2.6397 . 



7 Note that 7r/(2v / 2) = 1.1107 so that the bound is also fulfilled near x = 0. 
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5 Analytic structure in the complex plane 



5.1 The branch points 

The implicit equation (1.3) can be written in the form 

x 1 y + ix 

y = Arctan— = — Ln (5-1) 

y 2i y — ix 

where Arctan, Ln denote the general (multivalued) inverse tangent and logarithmic func- 
tions, respectively. One can immediately see that the function y = W t (x) is finite in any 
finite part of the complex x plane. Indeed, the only singularities of the expression on 
the r.h.s. of Eq. (5.1), as a function of x and y, are the logarithmic branch points at 
y = rkix. The latter cannot be satisfied simultaneously with Eq. (5.1) at any finite x. 
The derivatives of y = W t {x) have a rather simple structure (in implicit form): 

= dy = y = cos 2 i/ (5 2) 

dx y 2 + x 2 + x y + sinycosy 

,i d 2 y 2 x v 2y 3 

V = d^ 2 = ~(y 2 + x 2 + x) 2 ~ (y 2 + x 2 + xf 

„(") = *l = *j^y) ( 5 4 ) 

U dx n (y 2 + x 2 + x) 2n - 1 y ' ; 

where P n (x,y) is some polynomial in x and y. The first derivative diverges at the points 
where 

y 2 + x 2 + x = (5.5) 

together with the all higher derivatives. Because all the derivatives have the same poly- 
nomial (y 2 + x 2 + x) to an appropriate power in the denominator, they are either all finite 
or all divergent at the same point x. Therefore, in order to find all singularities of W t (x) 
it is sufficient to solve Eq. (5.5); this is equivalent to 

sin y cos y + y = (5.6) 

which involves only the function value y. One obvious solution is yo — 0; all other solutions 
of (5.6) are complex. Setting 2y = u + iv, where u and v are real, gives 

u -\- sin u cosh v = (5.7) 
v + cos u sinh v = 0. (5.8) 

We can restrict ourselves to u > (if y is a solution of (5.6) then — y is also a solution). 
Because cosh-y and smhv/v are always positive, we get the following constraints from 
Eqs. (5.7, 5.8) 

sin-u < 0, cosm < 0, (5.9) 
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which leads to 



(2n - 1) 7T < u < (2n - ir, n = 1,2, ... 
From Eqs. (5.7, 5.8) one gets an equation containing only the variable u: 

= \A 



tan if • arccosh 



anii 



' u 2 — sin 2 u 



(5.10) 



(5.11) 



which has exactly one solution u n in each interval n defined by Eq. (5.10). For large n 
one finds u n — b n — In (2b n ) /b n + . . . where b n = (2n — 1/2) tt and therefore 



±y n = \b n ± l -\n{2b n ) + 0\^Pi 

Xn - -\ M 2b n) -\±iK+o(^ 



(5.12) 
(5.13) 



which is in good agreement with the numbers in Table 2 obtained from the numerical 
solution of Eq. (5.11). 



n 








1 


-1.650611 ± 2.059981 i 


2.639705 


2.106196 ± 1.125364 i 


2 


-2.057845 ± 5.334708 i 


5.717853 


5.356269 ± 1.551574 i 


3 


-2.278470 ±8.522637i 


8.821948 


8.536682 ± 1.775544 i 


4 


-2.431122 ± 11.68877 i 


11.938917 


11.69918 ± 1.929404 i 


5 


-2.547991 ± 14.84580 i 


15.062869 


14.85406 ± 2.046852 i 


6 


-2.642706 ± 17.99809 i 


18.191069 


18.00493 ± 2.141891 % 



Table 2: The branch points x n of y = W t (x) and the corresponding function values y r , 



To investigate the type of the singularities we introduce new variables Ax and Ay in the 
vicinity of the singular point x n : 

x = x n + Ax , y = y n + Ay. (5.14) 

The implicit equation (1.3) can now be rewritten in the form 

Ax cos (y n + Ay) = (Ay cos Ay — sin Ay) sin y n + Ay sin Ay cos y n . (5.15) 
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If we go along a closed contour around the point y n in the complex y plane assuming that 
Ay is sufficiently small, then the corresponding path in the complex x plane will make a 
closed contour winding up twice as many times around the corresponding point x n . This 
can be demonstrated by inspecting the leading term in the expansion of Eq. (5.15) in 
powers of Ay: 

Axcosy n = {Ay) 2 cosy n + 0{{Ay) 3 ) (5.16) 
Therefore all singularities of W t (x) are of square root type : Ay = ±(Ax) 1//2 + ... or 

W t (x) ^ W t (x n ) ± Vx~=x^ . (5.17) 

Another way to obtain this result is to insert the following ansatz 

W t (x) ^ W t (x n ) +c(x- x n ) K (5.18) 

into the differential equation (3.5). Here k must be positive since we know that W t (x n ) 
exists. Assuming k < 1 we obtain in leading order 

Kc{x-x n ) K - x = ^-{x-x n )~ K (5.19) 

which determines k — 1/2 and c 2 = 1. 

5.2 The Riemann surface 

Since the function W t (x) has an infinite (countable) number of branch points, it has an 
infinite number of Riemann sheets as was already clear from Fig. 2. Therefore the question 
arises how to define these Riemann sheets in the most convenient way for any particular 
purpose. As we know the position and the type of all branch points, it only remains to 
define the cuts originating from these branch points. Given the cuts (we call any complete 
set of them a scheme), one then only needs to know how the different Riemann sheets are 
connected across the cuts. 

It turns out that the labeling of the Riemann sheets based on the inspection of the 
function on the real x-axis without a proper consideration of the branch points, as was 
done in Section 2, is not very illuminating when the whole complex rr-plane is concerned. 
Below we consider two other schemes of Riemann sheets which are better suited for the 
discussion of the singularities and analytical continuation across the cuts. In order to keep 
notations simple we prefer not to indicate the particular scheme explicitly in the notation 
w[ n \x); it will be defined in the context whenever the scheme is relevant. 

The first new scheme of Riemann sheets is demonstrated in Fig. 5 where the cuts 
originating from the branch points x n run to — oo along the lines of constant imaginary 
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sheet 


-1 


1 


2 


3 


n 


lim^+oo W t (n) (:r) 


. . . -tt/2 


tt/2 


3tt/2 


5tt/2 . . 


sgn(n)(|n| — l/2)n . . . 


lim^ +0 wl n) (x) 


-0 


+0 


7T 


2tt .. 


n{n — 1) 



Table 3: The scheme of Riemann sheets where the cuts extend to — oo as shown in Fig. 5. 




Figure 5: The Riemann sheets of y — W t (x) in the scheme where all cuts go to — oo: (a) 
sheet 1, (b) sheet 2. 

part Im x = Im x n , — oo < Re x < Re x n . Every Riemann sheet can be uniquely specified 
by the limit values at x — > +oo or x — > +0 as given in Table 3. 

There is a finite number of cuts on every sheet. The connection between different sheets is 
demonstrated in Fig. 5. On all sheets n except n = ±1, the function W t (x) is continuous 
and real on the real x axis. On the sheets n = ±1, the function is real for x > and 
purely imaginary for x < 0. 

Sending all cuts to infinity has, however, a drawback: it obscures the analytical prop- 
erties of W t (x) at x = oo. In fact, the points x = oo correspond to the regular points 
t — l/x — of the function W t (l/£). In the case when the analytical properties at infinity 
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Figure 6: The Riemann surface of y = W t {x): (a) sheet 1, (b) sheet 2, and (c) the images 
of different sheets in the complex y plane. The path a in plot (a) goes from sheet 1 to 
sheet —1 across the cut (Re xi, 0), the corresponding image of this path in the y plane is 
shown in plot (c). The paths j3 and 7 go from sheet 1 to sheet 2 across the cuts (Re x\,x\) 
and (xi,xl) respectively. 
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Figure 7: A modification of the main Riemann sheet of y = W t {x) shown in Fig. 6 (a): 
the cut between x\ and x\ is deformed so that it crosses the real x-axis by going over the 
branch point x — 0. 

are of interest, alternative schemes of the Riemann surface, such as shown in Fig. 6, have 
definite advantages. This second new scheme provides a convention of branches which 
is very suitable for the quantum mechanical example in Section 2, where the argument 
x was proportional to the inverse coupling constant so that the perturbative expansion 
became an expansion of W t around x = oo. In this scheme all cuts are finite in the com- 
plex x-plane: every branch point x n (n — 1, 2, . . . ) is connected with the corresponding 
complex conjugated point x* n , and every sheet, except for the sheets ±1, has a pair of such 
cuts. By deforming the cut connecting the branch points X\ and x* as shown in Fig. 7 
one can make it cross the real x-axis at x = 0— , so that the function is made continuous 
at both x > and x < 0, as in the case considered in Section 2. However, this comes at 
the price of having the branch point at x = superimposed on the cut between x\ and 
x*, which is not very convenient for the purpose of analytical continuation. 

The cut originating at x = occurs on every sheet in a finite interval: on the sheets 
±1 it goes along the real axis from zero to the cut connecting x\ and x\, on all other sheets 
it lies on the real axis between the two cuts connecting x n with x* n and x n _i with For 
every sheet n all the cuts are located within the limits \x\ < \x n \, Re x n < Re x < Re x n -±. 
The Riemann sheets in this scheme can be uniquely specified by the function values at 
infinity, which is a regular point on every sheet (see Table 4). 

This scheme will be used below in all references to different branches of W^^x). The 
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sheet 


-1 1 


2 


3 


n 


lim^oo W[ n) {x) 


. . . -tt/2 tt/2 


3tt/2 


5tt/2 . . 


sgn(n)(|n| — l/2)n . . . 



Table 4: The scheme of Riemann sheets with the finite cuts shown in Fig. 6. 



sheets n and — n differ only by the function sign: if y = W[ n (x), then — y = W\ {x) 
(see Eq. (2.19)). The connection between different sheets is shown in Fig. 6. 

The simplest way to investigate this connection is to begin with some path in the y 
plane (image) that connects two points on different sheets for the same x and to calculate 
the corresponding path in the x plane using Eq. (1.3). For any n = ±1, ±2, . . ., the sheets 
n and — n are connected to each other across the cut on the real axis between the points 
Re x n and Re x_„ (see the path a in Figs. 6 (a,c)). The cut between the branch points 
x n and x* n connects the sheet n with the sheet (n + 1) for n > or with the sheet (n — 1) 
for n < (the paths (3 and 7 in Figs. 6 (a,c) ). The structure of the cuts for all sheets 
with |n| > 2 is similar to that of sheet 2 (see Fig. 6 (b) ). 

The definition of W t (x) as an implicit function, Eq. (1.3), involves the function ytsmy 
which is single valued on the compactified complex plane. Therefore for any y there is one 
and only one value of x such that W t (x) = y. This allows us to determine the topology 
of any part of the Riemann surface by inspecting its projection onto the y plane. The 
topological structure of the Riemann surface of W t (x) turns out to be quite simple: any 
set of connected Riemann sheets with n — ±1, ±2, . . . , ±n max is homeomorphic to a two- 
dimensional disk as illustrated in Fig. 6 (c). For example, the part of the Riemann surface 
that consists of four sheets with n — ±1, ±2 is a manifold with a boundary homeomorphic 
to the one-dimensional sphere S 1 . The Euler characteristic 8 of this manifold is x — 1- 
This is equal to the Euler characteristic of a disk with a boundary. If this manifold is 
extended by adding pairs of sheets connected to its boundary (n = ±3, ±4, . . .) then 
the Euler characteristic remains unchanged. While much more could be said about the 
topological structure of the Riemann surface a further discussion is beyond the scope 
of this paper. We close by pointing out that, in a certain sense, the Riemann surface 
of W t (x) resembles an elaborate origami where a square of paper (which is topologically 

8 This is the topological invariant given hy \ = V — L-\-S where L, V, and S are the number of 
vertices, links and faces respectively in a triangulation of the manifold (sec e.g. Ref. [13]). For example, 
the triangulation with the vertices corresponding to the branch points and the links corresponding to the 
cuts in Figs. 6 (a,b) has V = 13, L = 16, S = 4. 
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equivalent to a two-dimensional disk) is folded 9 into a beautiful shape without paper cuts. 
The same property holds true for every function defined by an implicit equation involving 
any functions that are single-valued at finite x. 





Figure 8: (a) The imaginary part of the discontinuity of W^(x) across the cut — oo < 
x < 0, (b) the discontinuity of W[ 1 \x) across the cut connecting X\ and x\. 



5.3 The dispersion relation for W t (x) 



The analytic structure discussed above allows to derive the dispersion relation satisfied 
by W t (x) on the first sheet 



w t (1) w 

a 



7T 1 

2+^ 



dv 



AAv) 



At(v) 



U — Z 71 JO 

Re X\ , b — Im x\ . 



a + iv — z a — iv — z 



, (5.20) 
(5.21) 



Here A (u) and Ai(v) are determined by the discontinuities across the cuts on sheet 1 
(see Fig. 8): 



2i 
1 



A («) = - Wf'fw + iej-wr^-ie) = Im W^'iu + ie) , -oo < u < 0(5.22) 



Ai(u) = ^ I W t (1) (a + e + iu) 

9 The origami paper is not stretched, so the analogy ceases to be valid when we want to deform the 
disk into the Riemann surface under consideration. 



W t (1) (a-e + w) 



< v < b. 



(5.23) 
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6 Concluding remarks 



In the preceding Sections we have derived a number of results on the analytic structure, 
the numerical evaluation, series expansions and a simple quantum mechanical problem 
where the function W t {x) makes its appearance. While certainly not being exhaustive, 
these results demonstrate the amazingly rich structure and the fascinating properties of 
the W t -function which is defined by a simple, implicit equation. We hope that our selected 
findings may be useful for other researchers encountering this function. Of course, the 
methods used may be also adapted to related functions and indeed one might also consider, 
at the expense of convenience, more general implicit equations of which Eq. (1.3) is just 
a special case. For example, Wright has elaborated [14] on equations of the form 

{z + b)e z+a = z-b (6.24) 

and discussed the real solutions z to these equations as a function of a and b in Ref. [15]. 
Eq. (1.3) may be brought into this form, with a = and b = 2x, by writing the 
trigonometric function in its exponential form and setting W t = iz/2. Our results in 
Section 2 (and elsewhere) would then correspond to purely imaginary, special solutions 
of Wright's equation. 
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